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GRAPHITE  ABLATION  CHEMISTRY  NONEQUILIBRIUM  EFFECTS 


R.  L.  Baker 
The  Aerospace  Corporation 
Los  Angeles,  California 


Abstract 

The  implications  of  the  assumption  of  local 
solid-gas  phase  equilibrium  for  subliming  carbon 
species  for  graphite  ablation  calculations  in  an  air 
environment  is  investigated.  The  equilibrium 
assumption  is  eliminated  by  considering  the 
Knudsen-Langmuir  equation  at  the  interface  for 
each  carbon  specie.  Calculated  equilibrium  and 
nonequilibrium  results  are  compared  for  a very 
wide  range  of  flight  and  ground  test  environments. 
The  nonequilibrium  mass  addition  parameter  is 
always  less  than  the  equilibrium  value  and  the  non- 
equilibrium  wall  temperature  is  always  larger  for 
a given  environment.  Calculations  made  to  deter- 
mine the  convective  heat  flux  required  to  reach 
an  incipient  melt  temperature  of  3800°K  indicate 
that  the  required  flux  determined  from  an  equili- 
brium calculation  can  be  too  high  by  as  much  as 
200-300  percent  for  stagnation  enthalpies  less  than 
5000  Btu/lb.  Calculations  for  superorbital  reentry 
conditions  show  large  differences  in  the  mass  ad- 
dition parameter  B'  when  the  convective  heating 
rate  is  low  and  the  external  radiation  heating  level 
is  relatively  high.  Similar  large  differences  in  B' 
could  be  simulated  in  an  existing  ground  test 
facility  if  the  reported  external  radiation  heating 
level  could  be  increased.  Such  an  experiment 
would  provide  data  to  test  the  validity  of  present 
convective  heating  rate  blowing  corrections  in  a 
combined  heating  environment. 


Introduction 

In  the  formulation  of  problems  involving  in- 
terphase transfer  of  mass,  it  is  necessary  to  re- 
lato  the  concentrations  of  the  principle  mass 
transfer  species  in  one  phase  to  their  concentra- 
tions in  the  other  phase.  For  graphite  ablation 
calculations,  the  carbon  species  concentrations  in 
the  gas  phase  are  proportional  to  the  partial  pres- 
sures of  each  species  at  the  solid-gas  interface. 
Practically  all  graphite  ablation  calculations  are 
made  assuming  that  the  partial  pressure  p;  of 
carbon  species  i at  the  wall  is  equal  to  the  vapor 
pressure  p^  of  that  specie  at  the  wall  tempera- 
ture^-’). In  order  for  this  to  be  true,  the  solid 
and  the  gas  phase  must  be  in  equilibrium  with  one 
another  at  the  wall  temperature.  If  the  solid  and 
the  gas  phase  are  in  equilibrium  at  the  wall  tem- 
perature, then  the  same  numbor  of  molecules  are 
condonsing  from  the  gas  phase  into  the  solid  as 
those  vaporizing  from  the  solid  into  the  gas.  This 
means  there  is  no  net  transfer  of  mass  between 
phases,  therefore  we  know’  the  assumption  that 
p(  = p^  cannot  be  exact  for  finite  interphase  mass 
transfer  rates. 

The  purpose  of  this  work  is  to  formulate  the 
graphite  ablation  problem  without  making  the 
assumption  that  pj  = p'j'  and  to  determine  the  im- 
plications  of  this  simplifying  assumption  for  an  air 


environment  by  comparison  of  results  calculated 
using  both  methods.  In  the  remaining  discussion, 
the  problem  formulation  with  p,  assumed  equal  to 
p^  at  the  wall  temperature  is  simply  termed  equi- 
librium. When  this  assumption  is  not  made,  the 
term  nonequilibrium  formulation  is  used. 


Problem  Formulation 


The  Knudsen-Langmuir  Equation 


M) 


It  can  be  shown  from  kinetic  theory  ' that  the 
mass  of  gaseous  species  i striking  a unit  area  of 
wall  per  second  is  given  by 


0»; 


2itRT 


(1) 


_ is  the  wall  temperature,  OIL  is  the 


where  T^ 

molecular  weight  of  species  i and  p^  is  the  pres- 
sure. By  referring  to  Fig.  1,  we  find  expressions 
similar  to  that  above  can  bo  written  for  the  mass 
vaporizing  from  the  liquid  and  condensing  from  the 
gas  phase  per  unit  area  per  second.  Then  a simple 
interphase  mass  balance  gives  the  following  ex- 
pression for  the  net  transfer  of  mass  per  unit  area 
per  second  between  phases 


mi  = °i  y 2uRT  I 

v W 


(2) 


(4) 

This  equation  is  the  Knudsen-Langmuir  oquation  . 
The  pj'  and  pj  are  the  vapor  pressure  of  species  i 
at  the  wall  temperature  Tw,  and  the  partial  pres- 
sure of  species  i,  one  mean  free  path  from  the 
wall.  These  pressures  cannot  be  equal  to  one 
another  if  there  is  to  be  any  net  interphase  trans- 
fer of  mass.  The  coefficient  <»l  in  Eq.  (2)  is  the 
vaporization  coefficient.  From  kinetic  theory, 
the  maximum  mass  flux  vaporizing  per  unit  area 
of  wall  per  second  is  /2irR T^7  p V.  The  actual 
mass  flux  is  a ,/2rfRTu,  p^.  The  valuo  of  Oj 
must  be  determined  experimentally*'*).  In  ob- 
taining Eq.  (2),  it  has  been  assumed  that  the 
vaporisation  coefficient  and  tho  condensation  coef- 
ficient are  equal  (see  Fig.  1).  Equation  (2)  repre- 
sents a surface  mass  balance  on  a microscopic  or 
molecular  level. 

Mass  Conservation  Equation  - Binary  System 

The  principal  equations  used  m defining  tho 
present  problem  are  the  macroscopic  mass  and 
enorgy  conservation  equations.  As  shown  by 
Kubota^),  modeling  of  the  convective  and  diffu- 
sive mass  fluxes  at  tho  wall  (see  Fig.  2)  gives  tho 
following  relationship  between  the  mass  addition 
parameter  B1  and  the  mass  fraction  of  vaporizing 
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For  the  nonequilibrium  formulation,  the  pj  is  ob- 
tained from  the  Knudsen- Langmuir  relationship, 
Eq.  (2),  Thus,  if  Eqs.  (2)  and  (3)  are  combined, 
Pj  and  pj  are  related  by 
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Fig.  1.  Microscopic  Surface  Mass  Balance 
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where  p u C.T 
re  e H 


is  the  heat  transfer  coefficient. 
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For  m/o-iPj  « 1 , Eqs.  (1),  (2)  and  (4)  may  be 
combined  to  give  the  following  expression  for  the 


reduction  in  Bg^  due  to  nonequilibrium  effects 


eq 


(7a) 
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For  the  binary  system  presently  being  con- 
sidered, the  Wl  in  Eq.  (6)  can  be  written  in  terms 
of  0«i,  OW 2 » PI  and  pc.  This  expression  can  then 
be  combined  with  Eq.  (4)  to  give  the  following 
approximate  relationship  between  B'  and  B^  for 
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fig.  2.  Macroscopic  Surface  Mass  Balance 
species  at  the  wall,  Kt 


B'  = B’ 


eq 
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(7b) 


B1  = 


1 - K, 


(3) 


Kj  is  related  to  the  partial  pressure  at  the  wall 
PjWand  the  molecular  weight  0)1  j of  the  vaporizing 
specie  by  the  relationship  Kj  = pi*)Mj/p60>l  where 
pe  is  the  pressure  at  the  edge  of  tnc  boundary 
layer  and  (jjjis  the  average  aiolecular  weight  of 
the  gas  at  the  wall.  Making  this  substitution  into 
Eq.  (3)  and  simplifying  for  a binary  system  of 
components  1 and  2 results  in  the  following  equa- 
tion for  the  mass  addition  parameter  B': 


B'  -• 


' ^T 


-2--  1 


(4) 


In  the  equilibrium  formulation,  it  is  assumed  that 
Pj  - pv  leading  to  the  well-known  relationship 


From  Eqs.  (7a)  and  (7b)  wc  sec  that  the  nonequili- 
brium value  of  the  mass  addition  parameter  B1  is 
always  less  than  the  equilibrium  value  B^q.  The 
degree  of  departure  irom  equilibrium  n propor- 
tional to  m/pj,  or  PeueCpj/pe  and  inversely  pro- 
portional to  a j.  As  the  net  transfer  rate  of  mass 
between  phases  m increases  relative  to  the  rate 
of  mass  striking  the  wall  pj , nonequilibrium  et- 
fects  increase.  The  ratio  iVi/pj  is  proportional  to 
pcueCH/pe  so  that  maximum  nonequilibrium  ef- 
fects will  lie  seen  when  large  heat  transfer  coet- 
ficients  occur  at  low  pressures. 

Mass  Conservation  Equation  - Multiple 
Species 

The  above  equations  wero  restricted  to  a 
binary  system  with  no  surface  chemistry  other 
than  the  vaporisation  or  sublimation  process.  For 
graphite  ablation,  it  is  necessary  to  consider 
multiple  c rbon  species  at  the  wall  as  well  as 
chemical  roaction  of  the  carbon  with  the  free 
stream  gas.  Equation  (3)  for  graphite  ablation 
becomes^) 


For  unity  value  of  Lewis  and  Prandtl  numbers 


Equations  (7a)  and  (7b)  are  both  obtained  by 
expanding  in  a binomial  series  and  retaining 
the  first  two  terms. 


/A 


SKi  +fo  + fn 

B1  = 2 

1 - SK.  - F_. 
i N 
w 


where  Fq  and  F^  result  from  the  carbon -oxygen 
and  carbon -nitrogen  chemical  reactions,  respec- 
tively. The  mass  fractions  Kj  now  enter  as  a 
summation  over  all  carbon  species  i. 

The  functions  and  Fj^  are  given  by 


qc  “ qdiff  ' qconv  ‘ qr  + 


qc  " ^eue^H  " HW]  * ™[hw  ‘ ^*s]  * eaTw  + % 


qdiff  qc 


^e 


™C  "C2  C4 

fn  = [o«CN  kcn  kc2n2  + 0)!qj  ^ kc4n2 


where  Rn  for  air  is  0.232.  The  carbon-nitrogen 
compouncrmass  fractions  at  the  wall  were  present- 
ly calculated  assuming  gas  phase  chemical  equi- 
librium. * The  equilibrium  constants  were  ob- 
tained from  Dow  Chem^).  The  mass  fractions, 

Kj  , may  again  be  related  to  the  partial  pres- 
sures pj  and  the  molecular  weights  OWj.  So  that 
B1  may  be  written 


Sp/>VPe0»[Fo+FN) 
Pe<5Tl[i  - Fn]  - sp.oiij 


Equation  (8b)  is  applicable  in  general.  For 
tho  equilibrium  case,  the  pj’s  are  given  by  the 
individual  p^  at  the  wall  temperature.  For  the 
nonequilibrium  case,  pj  is  determined  from  Eq.  (6) 
written  for  each  of  the  i species.  It  should  bo_ 
noted  that  for  the  multiple  species  case,  the  Oil  in 
Eq.  (6)  depends  on  EpjOWj  so  that  an  iteration  pro- 
cedure is  required.  The  Oj's  used  in  Eq.  (6)  were 
the  "nominal  values"  givon  by  Dolton,  et  al,  i.e., 
Oj  = 0.  24,  <»2  = 0.  5,  Oj  = 0.  023,  <>4  = 0.  25  and 
05  = 0.0019.  All  the  calculations  presented  here 
were  obtained  using  the  carbon  species  thermo- 
chemical data  of  Palmer  and  SheleT*  to  calculate 
Pi  for  species  Cj  through  Cg  (see  Appendix  A). 

Energy  Conservation  Equation 

As  shown  in  Fig.  3,  tho  energy  conservation 
equation  at  the  interface  can  be  writton 

q = p u CU[H°-H  ] - m(H  -H  ] - <<rT4  + c 
Mc  e o Hl  w 1 w s'  w ^ext 


Fig.  3.  Macroscopic  Surface  Energy  Balance 


Nonequilibrium  effects  are  contained  implicit- 
ly in  the  energy  equation  through  effects  on  the 
ratio,  m/pgUgCj^,  i.e.,  B1  and  on  Hw  which  is 
evaluated  based  upon  p.  instead  of  on  pV. 

Primary  Dependent  and  Independent 
Variables 

Equations  (5)  or  (8)  and  (9)  provide  two  equa- 
tions in  the  primary  dependent  variables  B‘  and 
H^,.  Solving  these  equations  together  satisfies  the 
surface  mass  and  energy  conservation  equations 
simultaneously  and  provides  a unique  value  of  B‘ 
and  Hw.  For  the  equilibrium  case  both  B‘  and 
are  functions  of  the  edgo  pressure  p0  and  the  wall 
temperature  Tw.  If  the  steady  state  ablation 
assumption  is  made,  i.e.,  qc  = mfH^,-^],  Eq.  (9) 
may  be  written  _ 


where  the  substitution  PeueCj.j  H°  = qcv„  ) 

has  been  made.  The  term  H^is  the  bulk  ° 

enthalpy  of  the  solid  and  qcw  is  the  cold  wall  con- 
voctive  heat  flux. 


Usually  the  blowing  correction  Cjj/Ch0  ia 
written  as  a function  of  B’**.  Thus,  functionally 
Eq.  (10a)  may  bo  written 


F.  (p  , T ) F,(p  , T ) 

1 re  w 2 ro  w 

T 1 

H°  qcw 


where  q_  and  q8jrt  arc  the  net  conduction  heat  flux 
into  the  oody  ana  the  extornal  radiation  heat  flux; 
H°,  and  Hs  are  the  stagnation  enthalpy,  tho 
enthalpy  of  the  gas  at  the  wall  and  the  enthalpy  of 
the  solid  at  the  wall,  ih  is  the  total  mass  flux  of 
ablating  spocies  and  is  the  radiation  heat 

flux  away  from  tho  wall. 

$ 

This  simplifying  assumption  is  justified  by  the 
very  small  (<  5 percent  of  total)  surface  mass 
loss  associated  with  nitrogen  chemistry. 


Writing  the  onergy  balance  equation  in  this  way 
was  first  suggested  by  Kendall^8).  Equation  (10b) 
provides  valuable  insight  into  the  determination 
of  the  independent  variables.  The  edge  pressure 
pc  is  a primary  independent  variable.  The  '.vail 
temperature  Tw,  while  an  independent  variable,  is 

All  calculations  reported  in  this  work  were 
b^sed  upon  laminar  correction  to  the  Stanton 
number  to  account  for  surface  blowing  effects. 
The  relationship  used  is(^)  C H/Cu  = 
in(1.28B'  + D/1.28B1.  ^ 


3 


implicit;  the  other  variables  qcw  and  H°  are  not 
independent  but  rather  related  to  one  another  by 
Eq.  (10b).  If  the  value  of  B'  (and  therefore  of  T ) 
is  fixed  at  a given  pressure  pQ,  Eq,  (10b)  may  be 
written 


H 


A + rB_a  1 

,,o  q 


where  A and  B are  constants.  Thus,  the  selected 
value  of  B1  will  be  obtained  (at  the  pressure  pe) 
for  all  combinations  of  H°  and  qcw  satisfying 
Eq.  (10c).  A graphical  solution  procedure  based 
upon  this,  which  allows  the  determination  of  B1  for 
given  pe,  H°  and  q using  a single  graph  was 
given  by  Rindal,  et  al(°).  In  summation,  for  the 
equilibrium  case  B'  may  be  regarded  as  a func- 
tion of  pe,  H°  and  q . For  fixed  p , a given 
value  of  B1  is  obtained  for  all  H°  ana  qcw  satis- 
fying Eq.  (10c).  This  is  illustrated  in  Table  1. 


Table  1.  Equilibrium  B1  Values 


Pc 

latml 

*cw 

(Btu/ft2-scci 

H 

(Btu/lbl 

Tw 

<*K) 

B' 

0.1 

605 

100000 

3281 

0.420 

0.1 

647 

40000 

3281 

0.420 

0.1 

734 

20000 

i 

3281 

0.420 

0.1 

1000 

10000 

3281 

0.420 

0.1 

3646 

5000 

! 

! 3281 

0. 420 

0.1 

13112 

4400 

3281 

0. 420 

For  a pressure  of  0.  1 atmospheres,  B1  s 0.420 
for  all  of  the  H°,  qcvy  combinations  given  since 
thjy  all  satisfy  Eq.  (10c).  Note  that  the  wall 
temperature  is  tho  same  for  all  cases  also  as 
required.  Tne  .ndependont  variables  for  tho  non- 
oqu.1  rtum  case  are  discussed  in  the  next  section. 


Discussion  of  Results 
The  Nature  of  Ni  ncquilibnum  Effects 

To  initially  illustrate  differences  in  results 
calculated  using  the  equilibrium  and  the  nonequi- 
librium formulation,  it  is  informative  to  make 
such  a comparison  for  a senes  of  environmental 
parameters  for  which  tho  equilibrium  mass  addi- 
tion parameter  is  a constant.  If  the  previous  dis- 
cussion of  dependent  and  independent  variables  is 
recalled,  the  q , H°  conditions  in  Table  1 all 
give  a value  of  B'  (equilibrium  assumption)  of 
0.420.  Therefore,  let  us  compare  nonequilibrium 
calculated  results  wnth  the  equilibrium  results  m 


Table  1.  These  comparisons  are  shown  in  Fig.  4. 
The  calculated  equilibrium  results  are  all  given  by 
the  single  point  indicated  by  an  X on  the  equilibrium 
B'  curve.  Calculated  nonequilibrium  results  have 
a B'  less  than  the  equilibrium  value  and  a wall 
temperature  greater  than  when  solid-gas  equili- 
brium is  assumed.  The  reduction  in  B',  caused  by 
Pj  being  less  than  p^  [see  Eq.  (6)],  reaches  a 
reasonably  constant  value.  However,  the  surface 
temperature  continues  to  increase  over  its  equi- 
librium value  as  H°  is  decreased  and  q w in- 
creases, i.e.,  as  pgUgCrr  increases.  For  the 
equilibrium  case,  the  independent  variables  were 
shown  to  be  pe  and  all  combinations  of  fjcw  and  H° 
satisfying  Eq.  (10c).  From  Fig.  4 we  see  that  for 
the  nonequilibrium  case,  pe,  H°  and  qcw  are  all 
independent  variables  by  themselves. 

X ALL  EQUILIBRIUM  SOLUTIONS 
o qcw  = 605,  H*  = 100000 
C.  6 0 qcw  ^ 647.  H*  = 40000 

/ = 734>  H‘  = 20000 
0.5  / o 4CW  = 1°°°.  H*  = 10000 

Pe=0.1  atm  / . qcw  = 3646,  H*  = 5000 

0.4  A c qcw  = 13112,  H*  = 4400 


3100  3200  3300  3400  3500  3600  3700  3600 

TEMPERATURE  - *K 

Fig.  4.  Differences  Between  Equilibrium 
and  Nonequilibrium  Solutions 


Graphical  Solution  Procedure 

All  solutions  for  the  equilibrium  case  can  be 
represented  on  a single  graph  using  the  method  of 
Rindal,  et  al(8).  However,  since  there  are  three 
independent  variables  in  the  nonequilibrium  case, 
there  is  no  convenient  way  of  representing  all 
solutions  on  a single  graph.  For  a given  pres- 
sure pe  there  is  a unique  value  of  B*  and  Tw  (or 
specified  H°  and  1}  w as  show™  in  Fig.  4.  This 
suggests  that  for  given  pe,  a convenient  way  ot 
graphically  presenting  results  would  bo  to  plot 
lines  of  constant  H°  and  constant  f)  in  the  same 
B',  T plane  as  in  Fig.  4.  Such  a plot  for  an  edge 
pressure,  p,,,  of  one  atmosphere  is  shown  m 
Fig.  5.  For  this  pressure.  Fig.  5 gives  B'  values 
and  surface  temperature  fjr  q ranging  from  2000 
to  15000  3tu/ft2  -sec  and  for  the  11°  values  of 
5000  to  50000  Btu/lb.  Such  a plot  would  be  very 
useiul  for  instance  for  graphite  arc  jet  ablation 
results  at  constant  pressure  where  the  stagnation 
enthalpy  H°  and  the  cold  wall  heat  flux  q are 
varied  by  changing  the  arc  current  and  the  body 
geomotry. 
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Mass  Loss  and  Surface  Temperature  Differ- 
ences for  Typical  Reentry  Vehicla 

Environments  Fig.  6.  Nonequilibnum  Effect  on  Surface 

Recession  Rate 

Nonequilibrium  calculations  have  been  shown 
to  result  in  lower  values  of  B1  and  larger  wall 
temporaturcs.  The  surface  recession  rate  & is  * 
obtained  from  B'  through  the  relationship  *, 


where  is  the  bulk  density  of  the  ablating  mate- 
rial. Quantitative  results  showing  the  differences 
in  the  calculated  surface  recession  rate  5 and  the 
surface  temperature  Tw  for  a wide  range  of  re- 
entry environments  are  shown  in  Figs.  6 and  7, 

For  stagnation  enthalpies  of  2000,  5000  and  10000 
Btu/lb,  pressures  from  0.  1 to  100  atmospheres 
and  cold  wall  heating  values  from  300  to  75000 
Btu/ft2  -sec.  Fig.  6 shows  that  the  difference  m 
the  calculated  recession  rate  due  to  the  equili- 
brium assumption  is  nevor  greater  than  30  per- 
cent. Tho  corresponding  differences  in  tho  sur- 
face temperature,  shown  in  Fig.  7,  range  up  to  a 
maximum  of  about  600°K  obtained  when  tho  equi- 
librium surface  temperature  is  predicted  to  be 
about  3400°K,  whereas  the  calculated  nonoquili- 
bnum  temporature  is  about  4000°K.  In  both  Figs. 
6 and  7 we  see  that  the  largest  differences  occur 
at  lowor  pressures  (with  larger  PeUgCj^)  where  the 
ratio  of  the  interphase  mass  transfer  rate  to  the 
collision  frequency  at  the  wall  is  largest.  These 
differences  are  of  the  same  order  as  those  due  to 
the  assumption  of  alternate  carbon  species  ther- 
mochemistry models^),  Since  over  90  percent  of 
the  total  recession  occurs  at  pressures  greater 
than  10  atmospheres  for  typical  reentry  trajec- 
tories, it  can  be  seen  from  Fig.  6 that  the  differ- 
ence in  the  total  integrated  recession  will  bo  the 
order  of  5 to  10  percent  maximum.  Thus  the 
added  complexity  of  the  nonequilibnum  formula- 
tion is  not  justified  for  reentry  vehicle  nosetip 
surface  recession  predictions.  However,  a b00°h 
difference  in  the  surface  temperature  may  not  bo 
tolorable  for  thermostructural  design  purposes. 
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Fig,  7.  Nonequilibnum  Effect  on  Surface 
Temperature 


Graphite  Melt  Considerations 


The  question  of  whether  graphite  will  melt 
under  reentry  vehicle  convective  heating  environ- 
ments has  been  considered  by  Rindal  and  Powers^) 
and  by  Kratsch^®'.  These  authors  conclude  that 
for  high  beta  reentry  vehicle  trajectories  it  is 
possible,  for  small  nose  radii,  (1/4  to  1/2  in.)  to 
reach  an  assumed  melt  temperature  m the  range 
of  4000  to  4200°K. 


The  equations  described  in  this  work  are 
applicable  up  unt.l  the  time  that  tho  waii  tempera- 
ture reaches  the  melt  temperature.  When  the 
melt  temperature  is  reached,  additional  modeling 
of  the  melt  layer  on  the  surface  is  required. 
Equilibrium  and  nonequilibnum  calculations  were 
performed  to  determine  differences  m the  heat 
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flux  required  to  reach  an  incipient  melt  tempera- 
ture of  4200°K  at  a pressure  of  100  atmospheres 
for  stagnation  enthalpies  of  2000,  5000  and  10000 
Btu/lb,  For  H°  - 10000  Etu/lb,  the  required  heat 
flux  values  from  the  equilibrium  and  nonequili- 
brium  formulations  were  3000  and  2750  Btu/ft^- 
sec,  respectively.  For  H°  = 5000  Btu/lb,  the 
corresponding  required  heat  flux  values  were 
21700  and  6700  Btu/ft^-sec.  Thus  melting  could 
occur  at  this  stagnation  enthalpy  based  upon  the 
nonequilibrium  calculation.  However,  the 
required  equilibr.um  value  of  qcw  is  over  300 
percent  larger  and  could  not  be  obtained  at  this 
enthalpy. 

A melt  temperature  for  graphite  in  the  range 
4000-4200°K  is ,jiot  universally  accepted.  There 
is  some  old(H)  as  well  as  some  very  recent^' 
evidence  that  the  melt  temperature  is  about  3800°K 
and  that  the  triple  point  pressure  is  closer  to 
1 atmosphere  than  to  100  atmospheres.  There- 
fore, calculations  like  those  above  were  made  for 
an  incipient  melt  temperature  of  3800°K  for  pres- 
sures ranging  from  10  to  100  atmospheres.  The 
results  are  shown  in  Fig.  8.  The  difference  in 
the  calculated  heat  flux  required  to  attain  incipient 
surface  melting  is  quite  small  over  the  whole 
range  of  pressure  for  H°  - 10000  Btu/lb.  For 
H°  = 5000  Btu/lb,  the  equilibrium  prediction  for 
heat  flux  is  50  percent  too  high  at  lower  pres- 
sures. For  ground  test  environments  with  stag- 
nation enthalpies  in  the  range  2000-5000  Btu/lb, 
differences  of  over  200  percent  can  occur  even  at 
a pressure  of  100  atmospheres. 


Fig.  8.  Nonequilibrium  Effects  on  Melt 
Temperature  Environments 


1 am  indebted  to  Dr.  A.  G.  Whittaker  for  bring- 
ing this  reference  to  my  attention. 


These  results  indicate  very  large  differences 
in  the  predicted  heating  environment  required  to 
reach  incipient  melt  conditions  can  occur  when  the 
solid-gas  equilibrium  assumption  is  made.  Once 
the  melt  temperature  is  reached,  the  liquid  graph- 
ite has  a vaporization  rate  (B ’ ) much  greater  than 
the  solid  graphite^®*  H).  Because  of  this,  essen- 
tially 100  percent  of  the  liquid  melt  layer  is  vapor- 
ized and  very  little  is  stripped  off  as  liquid.  Thus, 
nearly  full  advantage  may  still  be  taken  of  the  very 
large  amount  of  heat  required  to  change  solid 
graphite  to  gaseous  species  in  the  boundary  layer. 
However,  when  the  surface  temperature  has 
reached  the  melt  temperature,  the  modeling  must 
be  changed  to  account  for  the  nresencc  of  liquid 
graphite  as  done  by  Kratsch-  No  attempt  was 
made  to  carry  out  this  type  of  modeling  in  the 
present  work. 

Radiation  Heating  Etfects 

Results  presented  thus  far  have  been  restricted 
to  a convective  heating  environment  for  ground  test 
or  flight  applications.  If  the  free  stream  velocity 
is  increased  from  reentry  vehicle  values  like 
25000  ft/sec  to  superorbital  velocities  like  50000 
ft/sec,  gas  cap  radiation  heating  to  the  solid  sur- 
face becomes  important  in  addition  to  the  convec- 
tive heating.  Graphite  ablation  predictions  for  an 
Apollo  type  vehicle  reentering  the  earth's  atmos- 
phere at  superorbital  velocities  have  been  made  by 
Bartlett,  et  al(^). 

Equilibrium  and  nonequilibrium  calculations 
were  made  for  the  range  of  convective  and  radia- 
tion heating  levels  givenU^l.  The  particular  con- 
ditions chosen  were  an  altitude  of  220000  ft  and  a 
velocity  of  45000  ft/sec  which  result  m a stagna- 
tion pressure  of  0.23  atmosphere  and  a stagnation 
enthalpy  of  40318  Btu/lb.  The  calculated  results 
are  summarized  in  Table  2.  The  surface  temper- 
ature difference  is  quite  small  over  the  whole 
range  of  cond.tions.  The  surface  recession  rate 
difference  is  small  until  the  condition  of  low  con- 
vective heating  rate  and  higher  radiation  heating 
rate  associated  with  larger  nose  radii  is  approached. 
For  these  heating  conditions  the  calculated  B^q  may 
be  nearly  twice  the  nonequilibrium  value,  with  both 
B'  and  Bg_  becoming  significantly  greater  than 
one,  i.e.,  massive  blowing.  Despite  the  large 
difference  in  B',  the  difference  m the  predicted 
recession  rate  s is  only  about  25  percent  because 
of  a compensating  effect  associated  with  the  blow- 
ing correction  to  the  convective  heating  rate. 

Thus,  there  appears  to  be  no  substantial  differ- 
ences in  either  the  wall  temperature  or  the  surface 
recession  rate  for  superorbital  reentry  conditions. 

The  above  statement  regarding  differences 
depends  intimately  upon  the  validity  of  the  blowing 
correction  to  the  convective  heat  trar.sler  rato  in  a 
radiation  augmented  heating  environment  where 
B'  is  boing  driven  to  very  large  values  by  the  ex- 
ternal radiation  heating.  This  type  of  coupling  was 
not  present  .n  the  data  upon  which  p.esent  blowing 
corrections  are  based.  It  would  thus  be  desirable 
to  obtain  new  data  in  which  B1  is  very  large  due  to 
external  radiation  heating.  In  Table  3,  calculated 


When  the  equilibrium  B'  is  much  larger  than  the 
nonequilibrium  B',  C^/Cf^  »»  much  less  for  the 
equilibrium  case.  Since  b - B'  qcw  C.,/C„  » 
the  two  effects  compensate.  ” 


b 


Tabic  2.  Superorbital  Reentry  Calculations 


(atm) 

H* 

(Btu/lb) 

(Btu/ft2-*ecl 

<Btu/u2_,ec) 

T -T 

w 

B' 

W 

B’ 

V* 

i 

0.23 

40318 

2503 

106 

157 

1.98 

0.232 

0.122 

0.23 

40318 

1823 

166 

125 

1.70 

0.225 

0. 124 

0.23 

40318 

1191 

290 

88 

1. 44 

0.220 

0.128 

0.23 

40318 

851 

426 

67 

1.41 

0.222 

0.130 

0.23 

40318 

596 

585 

51 

1.64 

0.236 

0.131 

0.23 

40318 

358 

825 

41 

3.32 

0.338 

0.146 

0.23 

40318 

237 

975 

42 

9.88 

0.746 

0.201 

Table  3.  Large  B1  in  a Ground  Test 
Environment 


(*tm) 

BtWib) 

*4U 

VT* 

<**> 

6* 

P’^-O* 

P‘ 

ki 

i 

•0.J2 

4698 

705 

1735 

ui 

0.86 

0.320 

0 20) 

*0.32 

4691 

705 

2193 

194 

1.18 

0,341 

0.203 

.0.32 

4691 

105 

3000 

213 

1.94 

0.389 

0.196 

•0.32 

4491 

705 

5000 

256 

5.84 

0.578 

0.195 

.0.32 

4491 

105 

6000 

276 

M.O 

0 809 

0 221 

'Conation*  I"  frwnt*  Of  WafcaH.w'* 
• KfMtM  4^,,  o*©*»  14  rtlutl 


B1  and  B^q  values  (or  several  radiation  heating 
levels  at  and  above  thoso  in  the  ground  test  ex- 
periments of  Wakefield^)  are  tabulated.  From 
this  we  see  that  questions  regarding  the  validity  of 
present  blowing  corrections  to  the  convective 
heating  rate  in  highly  coupled  i adiation-convective 
heating  environments  could  bo  explored  by  con- 
ducting now  ground  tests  with  radiation  levols  2 
to  3 times  highor  than  those  in  Wakefie».,,s  ex- 
periments. 


Summary  and  Conclusions 

Comparison  of  equilibrium  and  nonequilibrium 
calculated  results  for  the  same  environmental 
conditions  indicates  that  the  mass  addition  para- 
meter, B1,  ' s always  less  and  the  wall  tompera- 
turo  is  always  greater  for  the  nonoquilibrium  case. 
For  a given  pressure,  the  difference  in  B*  levels 
off  and  the  difference  in  wall  tomperature  con- 
tinues to  increase  as  increases  and  ^de- 
creases, i.o.,  as  the  heat  transfer  coefficient 
increases. 

For  the  nonequilibrium  case,  pe,  H°  and 
qcw  are  all  independent  variables.  A convenient 


graphical  solution  procedure  for  a given  pressure 
is  to  plot  lines  of  constant  stagnation  enthalpy  H° 
and  lines  of  constant  cold  wall  heat  flux  qcw  in  the 
B ’ , Tw  plane . " ' 

Comparison  of  calculated  equilibrium  and  non- 
equilibrium results  for  a wide  range  of  reentry 
vehicle  environmental  conditions  indicates  a maxi- 
mum difference  in  the  calculated  recession  rate  of 
30  percent  and  a maximum  difference  in  the  sur- 
face temperature  of  600°K.  The  largest  differ- 
ences occur  at  low  pressures  for  large  heat  trans- 
fer coefficients. 

If  the  graphite  melt  temperature  is  3800°K, 
then  nonequilibrium  calculations  should  be  used  to 
determine  the  heat  flux  required  for  incipient 
melting  for  stagnation  enthalpies  less  than  5000 
Btu/lb  since  differences  as  large  as  200  to  300 
percent  can  occur. 

Radiation  augmenting  of  the  convective  heating 
environment  can  cause  a dramatic  increase  in  B1. 
When  this  happens  differences  m B1  of  over  100 
percent  can  occur.  However,  the  corresponding 
difference  in  the  surface  recession  rate  is  only 
about  25  percent.  This  effect  occurs  for  super- 
orbital reentry  velocity  conditions  when  the  con- 
vective heating  rate  is  low  and  the  radiation  heating 
rate  is  high,  i.e.,  for  large  nose  radii.  Ground 
test  simulation  of  this  effect  is  possible  by  in- 
creasing the  external  radiation  level  of  tl.v  experi- 
ments of  Wakefield  2 to  3 times.  In  this  way, 
experimental  determination  of  the  validity  of  the 
blowing  correction,  Ch/Cj-jo,  could  be  made  for  a 
coupled  convective -radiation  heating  environment 
at  the  very  high  B1  values  of  interest. 


Appendix:  Thermochemistry  Data 

All  of  the  calculations  presented  here  were 
obtained  using  the  carbon  species  Cj  through  Cj 
thermochemical  data  of  Palmer  and  Shelofft).  The 
values  used  were  termed  in  that  vory  extensive 
survey  paper  the  "approximate  best. 11  Putting  the 
carbon  species  vapor  proesuro  equations  in  the 
form 


j v 

In  Pj  (atm)  = yjo^y  + B 


(A- 1 ) 


the  Palmer  and  Shclef  values  of  A and  B^'  are 
tabulated  in  Table  A-l. 


Table  A-l.  Constants  from  Palmer 
and  Shelof 


A B 


- 85715 

18.  69 

- 98363 

22.20 

- 93227 

23.93 

-150307 

31.  30 

-133087 

32.  71 
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